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PREFACE 


This  report  covers  work  accomplished  under  contract 
DAAG53-76-C-0134.  The  contract  was  for  the  survey  and  analysis 
accomplished  and  did  not  include  any  other  sub-tasks. 
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I.  INTRODUCTION 


I 


The  objective  of  this  effort  was  to  investigate  thermal  infrared 
(3-5  n and  8-14  1)  camouflage  measures  for  military  targets  in  foliated 
backgrounds.  Common  lightweight  materials  were  considered  as  the 
possible  emulators  of  the  foliated  backgrounds.  These  materials  were 
considered  to  be  cloth  like  material  with  very  low  thermal  inertia. 

Since  the  thermal  infrared  camouflage  problem  exists  throughout 
the  day,  the  infrared  signature  was  considered  for  the  emulator  and  the 
foliage  during  the  entire  day.  Both  analytical  and  experimental  evalua- 
tion of  the  foliage  signature  and  analytical  evaluation  of  a proposed 
emulator  were  examined. 


V 


Successful  evaluation  of  tft^thermal  infrared  signature  of  foliage 
for  a 24  hour  period  depends  upon  the  modeling  of  the  energy  exchange 
mechanism  of  the  foliage.  The  energy  exchange  of  foliage,  in  general 
terms,  cannot  be  modeled  successfully  because  of  the  many  independent 
variables,  however  the  energy  exchange  for  an  individual  element  of  the 
foliage,  such  as  a single  leaf,  may  be  modeled.  If  camouflage  materials 
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are  produced  which  emulate  the  infrared  signature  of  a 1 individual  leaf 
and  if  these  materials  are  properly  distributed  on  netting,  the  emulation 
of  background  foliage  should  be  possible.  For  this  reason,  the  major 
thrust  of  this  effort  was  toward  the  modeling  of  an  individual  leaf  rather 
than  a gross  foliage  background. 

Both  a leaf  and  the  leaf  emulator,  or  camouflage  material,  exist  in 
an  environment  which  exchanges  energy  with  them  by  two  basic  mechanisms, 
convection  and  radiation.  The  convection  mechanism  is  controlled  by  the 
external  flow  field,  wind  velocity  and  turbulence,  and  by  the  physical 
dimensions  of  the  object.  These  mechanisms  have  been  extensively 
investigated  and  are  reported  in  many  standard  texts  and  references, 
however,  the  results  for  the  specific  case  of  leaves  was  reported  by 
Parkhurst,  et  al.  in  reference  1, 

In  the  case  of  radiation  energy  exchange,  the  leaf  environment  is 
logically  divided  into  two  broad  radiant  energy  wavelength  regions.  These 
are  the  thermal,  wavelengths  of  about  2.  5 micrometers  to  20  micrometers, 
and  the  solar,  wavelength  of  about  0.  3 micrometers  to  2.  5 micrometers, 
regions.  Both  environments  have  been  extensively  investigated  and  may 
be  modeled  under  specific  meteorological  conditions.  However,  with  the 
very  large  number  of  meteorological  and  other  environmental  variables 
possible,  modeling  must  be  accomplished  with  average  conditions.  These 
average  conditions,  while  useful  for  comparison  of  IR  signatures  from 
camouflage  materials  and  leaves,  may  never  be  duplicated  by  natural 


II.  THERMAL  MODEL  FOR  LEAF 


Basic  Considerations 

The  transient  thermal  response  of  a leaf  is  governed  by  the  energy 
exchange  with  the  environment  at  the  leaf  boundaries  and  the  thermal 
capacitance  of  a leaf.  By  considering  a typical  leaf  thickness,  the  rate 
of  response  or  time  constant  is  of  the  order  of  300  seconds  [2].  Since 
the  time  periods  of  interest  in  this  study  are  much  longer,  the  leaf  was 
modeled  assuming  steady  state  occurs  in  each  time  period.  With  this 
basic  assumption  the  model  becomes  simplified  in  that  the  environmental 
input  may  be  equated  to  the  leaf  energy  loss  to  determine  the  leaf  temper- 
ature. 

Environmental  Conditions 

Convective  Environment:  Tho  convective  environment  for  a leaf 

consists  of  the  air  temperature  and  velocity.  Since  both  of  these  are 
meteorological  quantities  the  approach  taken  was  to  assume  a typical 
variation  for  diurnal  temperature  vaiiation  and  a constant  air  velocity. 
The  air  temperature  w?  s specified  over  the  twenty  four  hour  period  by 
the  ASHRAE  recommended  procedure  [3].  Maximum  and  minimum 
temperatures  were  chosen  for  the  t^  oe  season  of  interest. 

Wind  velocity  was  considered  constant  at  a value  input  into  the  pro- 


gram, Values  from  6.7  m/s  (15  rrph)  to  calm  were  examined. 


Radiation  Solar  Environment:  The  solar  energy  incident  upon  a leaf 
was  determined  using  the  ASHRAE  method  to  determine  the  intensity  of 
the  direct  solar  beam  and  the  diffuse  or  scattered  component  [3],  The 
solar  declination  for  any  day  was  used  as  a program  input.  From  this 
information  and  the  leaf  direction  parameters,  i.e.  the  angle  of  tilt  and 
the  azimuth  angle,  the  solar  energy  incident  on  the  leaf  upper  surface  was 
evaluated.  Reflected  solar  radiation  was  assumed  to  be  negligible. 

Longwave  Radiation  Environment:  Longwave  radiant  energy  impinges 
on  the  leaf  from  three  sources.  First,  from  the  warm  humid  atmosphere; 
second  from  the  ground  below  the  leaf;  third  from  any  object  in  view  of 
the  leaf.  The  radiant  energy  originating  in  the  atmosphere  is  a function 
of  the  atmospheric  dry  bulb  temperature  and  the  moisture  content.  This 
radiant  flux  is  approximated  by  several  authors  using  an  equation  of  the 
form:  [4] 

q1'  m = a T^A  + BVe*)  1 

atm  a 

in  which 

2 

q"  t is  the  radiant  flux  in  cal.  /cm  /min. 
atm 

T&  is  the  absolute  air  temperature  in  °K 

e is  vapor  pressure  of  the  water  in  the  air  in  millibars 

a is  the  Stefan  Beltzmann  Radiation  constant 

A and  B are  constants  varying  from  author  to  author  in  \ approximate 

ranges  of  0.4  £ A $ 0.  75,  0,  047  s B ^ 0.  08 

The  long  wavelength  energy  flux  on  the  leaf  from  the  ground  is  given 
by:  i 
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q"  = « a T F , 

gr  gr  gr  t - g 2 

in  which, 

2 

q"  is  the  radiant  flux  in  cal.  /cm  /min. 

gr 

T is  the  absolute  temperature  of  the  ground, 
gr 

c is  the  long  wavelength  emittance  of  the  ground  i.  o,  , * (T  ), 
gr  gr  gr 

F is  the  fraction  of  the  ground  seen  by  the  leaf,  i.  e.  one  for  the 
i -g 

bottom  of  a horizontal  leaf,  £ for  one  side  of  a vertical  leaf,  etc. 

The  third  component  of  the  long  wavelength  radiant  energy  is  a variable 


dependent  upon  the  individual  case  studied.  For  the  purpose  of  this  work, 
this  term  was  assumed  to  be  negligible.  This  assumption  was  made  on 


the  basis  that  an  emulator  would  react  to  surrounds  exactly  like  the  leaf 
if  the  emulator  reacts  to  the  simplified  surrounds  like  the  leaf. 

In  order  to  evaluate  radiant  input  to  the  leaf  from  the  ground,  the 
ground  temperature  must  be  determined.  Modeling  of  the  ground  to  obtain 
this  temperature  requires  a transient  model.  In  this  case  the  ground  was 
modeled  as  a semi-in  inite  slab  with  constant  thermal  properties.  Properties 
used  were  obtained  from  reference  5.  One  of  the  most  uncertain  variables 
in  this  sub-analysis  was  the  convective  heat  transfer  from  the  ground 
surface  to  the  air.  This  was  finally  obtained  from  reference  6 using 

reasonable  vegetation  heights.  The  expression  used  was 
# 

h = (Cf/2)pCp(Pr)  u 3 

Budyko  [7]  presents  results  of  a study  of  drag  caused  by  vegetation  on  the 


earth's  surface  and  proposes  a drag  coefficient  as  follows: 


C /2  = [2.51n(y/Z  ) + 5]' 
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in  which: 

2 

h is  the  convective  coefficient  in  cal.  /cm  /min.  - *C 

Cf  is  the  drag  coefficient  depending  upon  surface  roughness, 

3 

p is  the  density  of  air  g/cm  , 

C is  the  specific  heat  of  air  cal/gm-  UC, 

P 

P is  the  Prandtl  number  of  air, 
r 

U is  the  wind  velocity  cm/min, 

y is  the  height  of  wind  velocity  measurement  assumed  to  be  76  cm, 

Zq  is  the  height  of  vegetation  in  cm. 

Leaf  Energy  Balance: 

Radiant  Input:  The  radiant  energy  input  to  a leaf  was  divided  into  the 
two  parts,  solar  energy  input  and  long  wavelength  energy  input.  The  solar 
normal  or  beam  intensity  was  obtained  from:  [3] 

q"SN  = A/exp(f3/ sin  a)  5 

in  which 

q"  is  the  normal  solar  intensity  v>/m^ 

O IN 

a is  the  solar  altitude  in  degrees 
A is  the  apparent  solar  irradiation  at  air  mass^O 
3 is  the  atmospheric  extinction  coefficient 


The  actual  energy  incident  on  the  leaf  consists  of  the  solar  beam  radiation 

projected  on  the  leaf  surface  plus  the  scattered  solar  energy.  Scattered 

solar  flux  was  determined  from:  [3] 

q1'  = Cq"  F 6 

q SD  4 SN  l-s 

in  which 
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q"  is  the  diffuse  or  scattered  solar  irradiation, 

Di-J 

q"g^  is  defined  in  equation  5 

C is  a constant  which  depends  on  atmospheric  dust,  moisture,  and 
air  mass, 

F is  the  fraction  of  the  sky  seen  by  the  leaf,  configuration  factor 
X-  s 

from  the  leaf  to  the  sky. 


These  two  energy  fluxes  were  assumed  to  be  incident  on  the  leaf  upper 

surface  and  the  energy  input  was  taken  to  be  their  sum  multiplied  by  the 

solar  absorptance  of  the  leaf  upper  surface. 

q"  = (q"  + q"  ) a 7 

H TOT  SD  h SN  stu 

Values  for  the  solar  absorptance  of  leaves  are  available  in  many  sources 

but  the  values  used  were  from  Birkebak  and  Birkebak.  [8 ] 

Long  wavelength  energy  input  to  the  leaf  from  the  air  and  the 

ground  (qM  ) as  given  in  equations  1 and  2 were  treated  as  follows.  The 

energy  input  from  the  atmosphere  was  assumed  to  be  the  same  for  both 

the  upper  and  lower  surfaces  of  the  leaf  and  the  ground  energy  input  was 

assumed  to  be  to  the  leaf  lower  surface  only.  Thus  the  total  energy  input 

per  unit  area  of  leaf  was  a q"  , + cr  q"  . +Qf,  q"  where  a,  is  the 

r tu^  atm  t£  atm  tjj^  gr  tu 

long  wavelength  absorptance  of  the  leaf  upper  surface  and  a is  the  long 

ti 

wavelength  absorptance  of  the  leaf  lower  surface. 

Convective  Energy  Input:  The  convective  energy  input  was  evaluated 
from  the  standard  convective  heat  transfer  expression 


q"  = h (T  -T.  ) 

conv  c atm  leaf 
in  which 
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q"  is  the  convective  heat  flux  in  cal.  /min.  -cm 

conv 

2 

h is  the  convective  heat  transfer  coefficient  in  cal/min-cm  -°C 
c 

T is  the  atmospheric  temperature  in  “C 
atm 

T,  , is  the  leaf  temperature  in  °C 
leaf 

Values  of  h^  were  obtained  using  the  procedures  of  reference  1.  These  are 
based  on  standard  equations  for  free  and  fcrced  convection  and  therefore 
require  the  evaluation  of  the  following  dimensionless  parameters: 

(1)  The  average  Nusselt  Number 

Nb  = h L/k.  9 

c 

(2)  The  Reynolds  Number 

R j = UpL/n  10 

J-1 

(3)  The  Grashof  Number 

GrL  = Pgp2L  At/p,2:  11 

in  which 

L is  the  effective  leaf  dimension  in  cm 
U is  the  wind  velocity  in  cm/min 

2 

g is  the  acceleration  of  gravity  in  cm/min*" 
k is  the  thermal  conductivity  of  air  in  cal.  /cm-min-“C 
3 is  the  temperature  coefficient  of  volume  expansion  in  cm 
la,  is  the  absolute  viscosity  of  air  in  gm/cm-min 

3 

p is  the  density  of  air  in  gm/cm 

The  standard  correlations  using  the  dimensionless  parameters  are 


(1)  For  free  convection  from  vertical  plates 
Nu  = 0.480  Gr  * 
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(2)  For  free  convection  from  the  upper  surface  of  a horizontal  plate 
warmer  than  air  or  to  the  lower  surface  of  such  a plate  cooler 
than  air: 

__  1 

Nu  = 0.497  Gr  4 13 

.Li 

(3)  For  free  convection  from  the  lower  surface  of  a warmer  than  air 

horizontal  plate  or  to  the  upper  surface  of  a cooler  than  air 

horizontal  plate: 

— 1 

Nu  = 0.  249  Gr  4 14 

Li 

(4)  For  forced  convection  to  or  from  a plate  having  a uniform 
temperature: 

NU  = 0.  595  Re.  ^ 15 
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Fig.  1.  Values  of  R for  ten  real  leaves.  X = flow  perpendicular  to 
stem;  L = flow  parallel  to  stem;  FORC.  = forced  convection  in 
laminar  flow;  FREE  = free  convection  in  laminar  flow;  TURB.  = 
forced  convection  in  turbulent  flow.  (Reproduced  from  Ref.  1.  ) 
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A decision  as  to  whether  the  flow  field  was  strong  enough  to  cause  forced 

Z 

convection  was  made  by  comparing  the  values  of  Gr  and  (Re  ) . If 

Li  L 

2 

(Re  ) was  larger  than  Gr  the  forced  convection  equation  was  used. 

Values  of  L for  use  in  these  expressions  were  obtained  from  L = RL 

max 

in  which  R was  obtained  from  figure  1 and  L was  the  average  maximum 

max  * 

dimension  of  the  leaf  under  consideration. 

Transpiration  Energy  Loss:  The  energy  loss  from  a leaf  by  transpira- 
tion of  the  leaf  moisture  was  extensively  studied  by  several  authors.  In 
this  work,  the  results  and  methods  presented  by  Gates  in  references  9 and 
10  were  used.  The  expression  for  transpiration  energy  loss  is 


r + r 
l a 


fg 
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in  which 

2 

q"e  is  the  transpiration  energy  flux  in  cal.  /cm  -min. 

pg(T^)  is  the  density  of  saturated  water  vapor  at  the  leaf  temperature  in 
gm/cm  , 

p (T  ) is  the  density  of  saturated  water  vapor  at  the  air  temperature  in 
g a 

gm/cm  , 

h is  the  latent  heat  of  vaporization  of  water  at  the  leaf  temperature 

cp  is  the  relative  humidity  of  the  air, 

r is  the  internal  leaf  diffusion  resistance  in  min/cm 
l 

r is  the  boundary  layer  resistance  given  by  equation  17. 


r = k, 
a 2 


w°-20l0*35 


u 


0.  5S 
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in  which 

W is  the  leaf  dimension  transverse  to  the  wind  in  cm, 

I-i  is  the  leaf  dimension  in  the  direction  of  the  wind  in  cm, 

U is  the  wind  velocity  in  cm /min 

1.  55 

is  the  dimensional  constant  of  0.247  cm/min. 

Values  for  the  internal  diffusion  leaf  resistance  are  presented  by  Gates 
for  several  common  leaves.  These  values  are  presented  as  constants 
although  it  is  known  that  water  stress  or  high  environmental  temperatures 
cause  these  resistances  to  change.  In  this  work  constant  leaf  resistance 
values  were  assumed. 

Leaf  Radiant  Energy  Loss:  The  radiant  energy  loss  of  the  leaf  was 
calculated  assuming  the  upper  and  lower  surfaces  had  the  same  emittance. 
This  results  in 

q"  =2®aT4  18 

r l 

in  which 

q"^  is  the  radiant  loss  per  unit  area  of  leaf  in  cal.  /cm^-min 

e is  the  emittance  of  the  leaf  at  the  leaf  temperature 

T is  the  absolute  leaf  temperature  in  °K 
l 

The  leaf  was  assumed  to  have  negligible  thermal  mass,  therefore,  the 
sum  of  the  energy  gain  was  set  equal  to  the  energy  loss  to  calculate  the 


leaf  temperature. 


Leaf  Radiance  in  the  3-5  and  8-14 
Micrometer  Wavelength  Region 

3-5  Micrometer  Wavelength  Region:  The  radiant  flux  from  a leaf  in 
the  3-5  micrometer  wavelength  range  was  assumed  to  consist  of  three 
components.  These  components  were;  (1)  energy  emitted  as  a function  of 
the  leaf  temperature;  (2)  energy  reflected  off  the  leaf  upper  surface  from 
incident  solar  energy  and,  (3)  solar  energy  reflected  off  of  the  ground 
and  transmitted  through  the  leaf.  The  components  were  calculated  as 


follows: 


R (Thermal)  = so  F T 
3-5  1 i. 

where 


2 


R^  (.(Thermal)  is  the  first  component  cal.  /cm  /min. 


c;  is  the  leaf  emittance 


T is  the  leaf  temperature,  °K 


Fjis  the  fraction  of  energy  radiated  between  3 and  5 micrometers  by  a 


Plankian  radiator. 


R,  _(Sol?.r)  = p F_G  cos0  20 

3-5  su  2 s 

where 

2 

R,  (Solar)  is  the  second  component,  cal.  /cm  /min. 

3-D 

p is  the  solar  reflectance  of  the  leaf  upper  surface 
su  rr 

F_  is  the  fraction  of  the  solar  insolation  which  is  in  the  3-5  micrometer 
2 

wavelength  range 

2 

G is  the  solar  insolation  in  the  3-5  micrometer  range,  cal.  /cm  /min. 
s 

» 

8 is  the  zenith  angle 


R,  -(Solar  Reflected)  = t p__  F_G  cosQ  21 

J-S  i aCi  c s 

where 

2 

R-  -(Solar  Reflected)  is  the  third  component,  cal.  /cm  /min. 

3-5 

t . is  the  leaf  transmittance, 
j l 

The  sum  of  these  three  would  be  the  energy  flux  at  the  leaf  surface  which 
would  be  detected  by  a 3-5  micrometer  wavelength  energy  sensor  system. 

8-14  Micrometer  Wavelength  Region:  The  radiant  flux  from  the  leaf 
in  the  8-14  micrometer  wavelength  range  was  assumed  to  be  from  the 
same  three  components  as  used  in  the  3-5  micrometer  wavelength  range. 
In  calculating  the  values,  the  only  change  in  equations  19,  20,  and  21  are 


the  fractions  and  F^.  Thus  the  same  equations  were  used  with  different 
values  for  these  fractions. 


III.  EXPERIMENTAL  VERIFICATION  OF  THERMAL  MODEL 
Experimental  Procedure:  In  order  to  test  the  validity  of  the  thermal 
model  used,  an  experiment  was  designed  and  run.  This  consisted  of 
instrumenting  leaves  on  two  living  trees,  measuring  the  leaf  temperature 
and  the  environmental  conditions  for  the  leaves. 

Leaf  Measurements:  A nearly  horizontal  leaf  on  a common  burr  oak 
(Quercus  Macrocarpa),  a vertical  leaf  on  a silver  maple  (Acer  Saccharnum) 
and  a horizontal  leaf  on  the  silver  maple  were  instrumented.  In  order  to 
minimize  the  effect  of  the  measurement  probes,  small  gage  (40  gage) 


thermocouples  were  installed  on  the  lower  surface  of  the  leaves.  These 
thermocouples  were  shielded  from  direct  solar  insolation  by  being  on  the 
lower  leaf  surface  and  were  in  contact  with  the  leaf  for  about  2 cm  which 


14 


should  reduce  thermocouple  conduction  effects.  Temperatures  were  con- 
tinuously recorded  on  a 12  point  Leeds  and  Northrup  Speedomax  W recorder. 
As  a check  on  the  measured  temperatures,  total  radiation  pyrometric 
temperatures  were  measured  with  a Barnes  PRT-10  radiometer. 

Ground  Temperature  Measurement:  Ground  temperatures  at  depths 
of  5 and  6.4  cm  were  made  using  24  gage  thermocouples.  Surface 
temperature  could  not  be  measured  directly,  due  to  short  grass  on  the 
surface,  therefore,  the  radiometric  temperature  was  measured  with  the 
Barnes  PRT-10  radiometer.  Actual  ground  temperatures  were  estimated 
from  the  radiometric  measurement  by  assuming  a ground  emittance 
consistent  with  the  actual  ground  surface. 

Environmental  Conditions:  The  environmental  variables  needed  for 
this  study  were  dry  bulb  air  temperature,  wet  bulb  air  temperature, 
wind  velocity  and  solar  insolation.  Dry  bulb  temperatures  were  measured 
continuously  using  a shielded  thermocouple  and  wet  bulb  temperature  was 
measured  each  30  minutes  using  a hand  sling  psychrometer.  Wind  velocity 
was  measured  at  1.  1 meters  above  the  ground  using  a ball  and  cup  anemo- 
meter. Solar  insolation  was  measured  continuously  using  an  Eppley  8-48 
black  and  white  pyranometer.  Experimental  data  is  summarized  in  Table  I. 
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Results  of  Experimental  Verification 

The  thermal  model  was  used  to  calculate  the  leaf  temperatures  and 
ground  temperature  for  the  72  hour  experimental  period.  Measured  values 
of  air  temperature,  relative  humidity,  wind  speed,  and  solar  insolation 
on  a horizontal  plane  were  used  at  one  half  hour  intervals  to  calculate 
leaf  and  ground  temperatures.  The  calculated  values  of  leaf  and  ground 
temperatures  and  the  experimental  values  are  shown  in  figures  2 through 
7.  Each  figure  is  for  one  full  days  data,  i.  e.  , for  the  entire  diurnal 
cycle.  The  figures  are  in  pairs,  that  is  figure  2 and  3 are  for  4 September 
1976,  figure  4 and  5 for  5 September  1976,  figures  6 and  7 for  6 September 
1976.  The  leaf  temperature  measured  and  calculated  are  shown  on  figures 
2,  4,  and  6.  Calculated  and  measured  ground  temperatures  along  with 
measured  air  temperatures  are  shown  in  figures  3,  5,  and  7, 

In  order  to  assess  the  accuracy  of  the  thermal  model,  the  errors  were 
statistically  examined.  This  analysis  indicated  the  mean  error,  i.  e.  , the 
measured  leaf  temperature  minus  the  calculated  leaf  temperature  to  be 
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0.  21  degrees  Celsius  with  a standard  deviation  of  the  errors  of  4.  00 
degrees.  From  this  information  the  leaf  temperature  calculated  has  a 
mean  error  with  95%  confidence  of  -0.  46  or  +0.  89°  Celsius. 

IV ^ CAMOUFLAGE  MATERIAL  REQUIREMENT 
Leaf  Temperature  Emulator:  The  thermal  model  prepared  was  used 

I 

to  determine  the  necessary  properties  of  a material  which  would  emulate 
leaves  thermal  response  during  one  day.  A study  of  the  leaf  response  and 
characteristics  indicated  that  ^n  emulator  could  be  produced  if  the  material 
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emittance  and  reflectance  could  be  controlled.  As  an  example  of  this 
figure  8 shows  the  thermal  response,  as  predicted  from  the  thermal 
model,  for  a maple  leaf,  an  aspen  leaf,  and  a material  designed  to 
emulate  the  maple  leaf.  This  material  would  be  a material  with  low 
thermal  mass,  clothlike,  with  solar  absorptance  of  0.47  and  long  wave- 
length emittance  of  0.  95.  As  can  be  seen  from  figure  7,  this  material 
would  satisfactorily  emulate  the  diurnal  temperatures  of  maple  leaves 
but  would  not  emulate  the  diurnal  temperatures  of  aspen  leaves.  This 
is  because  the  transpiration  rate  for  an  aspen  leaf  is  greater  than  the 
transpiration  rate  of  a maple  leaf.  An  emulator  for  the  aspen  would 
have  different  properties  which  could  be  found  using  the  thermal  model 
presented. 

Camouflage  Material  Detectability:  In  order  to  assess  the  detect- 
ability of  the  leaf  emulator  with  a maple  leaf  background,  the  total  radiant 
energy  leaving  the  emulator  and  the  leaf  were  evaluated  for  the  3-5 
micrometer  and  8-14  micrometer  wavelength  bands.  These  values  are 
plotted  in  figures  9 and  10.  Within  the  accuracy  of  the  model,  the  leaf 
radiance  and  the  emulator  material  are  identical.  This  indicates  that  the 
camouflage  role  of  the  material  would  be  well  fulfilled  in  the  3-5  and  8-14 
micrometer  wavelength  ranges. 

Visual  Camouflage  Problem:  The  solar  absorptance  of  the  camouflage 
material  which  was  used  to  obtain  the  results  shown  was  0.47.  The  solar 
absorptance  of  a material  is  a function  of  the  spectral  reflectance  over 
the  wavelength  range  in  which  the  suns  energy  reaches  the  earths  surface. 
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This  wavelength  range  is  approximately  from  0.  3 to  2.  5 micrometers. 
Since  the  visual  range  is  overlapped  by  the  solar  range,  the  visual 
reflectance  is  not  independent  of  the  solar  reflectance.  This  causes  a 
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camouflage  material  suitable  for  a leaf  emulator  to  be  more  reflective 
than  the  typical  visual  camouflage  material.  [11]  The  ideal  material 
to  be  used  for  both  visual  and  thermal  emulation  of  leaves  would  be  one 
which  had  a visual  reflectance  around  0.  3 or  lower  and  a solar  reflectance 
of  0.  53.  Such  materials  were  reported  in  reference  11  for  the  near 
infrared,  i.  e.,  0.7  to  1.2  micrometers.  Values  reported  are  shown  in 
Table  II. 


Table  II.  Infrared  Reflectance  for  Camouflage  Cloth  (Reference  11) 

(From  National  Military  Establishment 
Specification  JAN-C-765) 


Fabric 

Color 

(No.) 


Color 


Infrared  Reflectance  Percentages 
Relative  to  Magnesium  Oxide 
(Minimum  %) (Maximum  7o) 


1 

Light  green 

37.  0 

57.  0 

2 

Dark  green 

37.  0 

57.  0 

; 3 

Sand 

24.  5 

100 

? 4 

Field  drab 

24.  5 

57.  0 

* 5 

Earth  brown 

24.  5 

57.  0 

t 6 

Earth  yellow 

24.  5 

100 

i 8 

Earth  red 

24.  5 

57.  0 

9 

Olive  drab 

24.  5 

57.  0 

10 

Black 

0 

24.  5 

11 

White 

57.  0 

100 

12 

Forest  green 

24.  5 

100 

‘OYV\;. . 

13 

Desert  sand 

24.  5 

100 

However,  it  was  reported  that  the  production  of  such  materials  was  quite 


difficult.  Similar  work  was  reported  in  reference  12 
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Transparent  Cover  Analysis:  Since  it  may  be  difficult  to  obtain  the 
desired  visual  and  thermal  infrared  camouflage  with  a single  material, 
the  possibility  of  using  a visually  transparent  cover  over  an  opaque  material 
was  examined. 

Solar  Radiation,  Sky  Radiation 
Atmospheric  Radiation 


-f- 


Convection  Transfer 


/ Ground  Radiation 

Basis  for  Transparent  Cover  Analysis 
Figure  11 


^-Transparent  Sheet 
^ -Opaque  Sheet 


/ 


The  basic  system  analyzed  is  shown  in  figure  11.  A transparent  sheet  not 
in  contact  with  an  opaque  sheet  with  radiant  and  convective  heat  transfer 
was  considered.  The  energy  exchange  between  the  sheets  and  the  atmos- 
phere was  analyzed  by  consideri-.g  the  convective,  radiation  in  solar 
wavelengths  and  radiation  at  long  wavelengths  energies  as  uncoupled 
variables.  The  basic  radiant  energy  quantities,  considered  as  shown  dia- 
g ■ aromatically  in  figure  12.  In  this  figure  the  fluxes  indicated 


A 


\ j She*it  1,  Transparent 

1 

1 J 
1 S 2 

i 

T q" 

9 . SOl 

1 

i J 

s2 

' t2 

J J 

5i  i h 

4 

^ Sheet  2,  Opaque 

*2Eb. 


| tM  ground 


Radiant  Fluxes  for  Transparent  Sheet  Analysis 

Figure  12 
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Using  these  values,  energy  balances  on  sheet  1 and  2 including  convective 

energy  transfer  to  the  surrounds  and  between  the  two  sheets  results  in 

two  coupled  non-linear  equations.  These  equations  are  functions  of  the 

environmental  parameters;  q"go}>  ambient  air  temperature,  ground 

temperature,  and  the  radiative  properties  of  the  two  sheets.  In  order  to 

assess  the  possible  usefulness  of  the  transparent  outer  sheet  an  inner 

sheet  with  solar  absorptance  of  0.6  (dark  green)  was  considered.  Using 

ambient  air  Temperature  of  30  “C,  relative  humidity  of  50%,  ground 

2 

temperature  40  °C,  q"  of  800  w/m  , wind  speed  of  2.5  miles  Der  hour 

sol 

the  dark  green  material  temperature  calculated  was  41  °C.  With  a trans- 
parent sheet  over  the  dark  green  material  the  calculated  temperatures  of  the 
transparent  sheet  (T^)  and  the  dark  green  sheet  (T^)  with  several  different 
transparent  sheet  properties,  is  given  in  Table  III.  Notice  the  transparent 


Table  HI 

Sheet  Temperatures  Calculated 


Sheet 

Reflectance 

Sheet 

Transmittance 

Sheet 

Absorptance 

Tj. 

•c 

°C 

0.  10 

0.  80 

0.  10 

37 

46 

0.  05 

0.  80 

0.  15 

37 

46 

0.  05 

0.  75 

0.  20 

39 

46 

0.  10 

0.  70 

0.  20 

40 

45 

sheet  temperature  runs  from  37  to  40°  Celsius  where  the  uncovered 
material  temperature  was  41  “Celsius.  This  indicates  the  transparent 
sheet  over  the  camouflage  material  might  be  useful  if  the  proper  material 
cannot  be  obtained.  Limited  experimental  results  for  this  type  system  are 
reported  in  Appendix  B. 
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V.  CONCLUSIONS 

The  thermal  model  prepared  to  emulate  foliage  satisfactorily  predicts 
the  diurnal  temperature  of  leaves.  This  model,  like  all  models,  finally 
depends  on  the  quality  of  the  boundary  or  driving  functions.  In  this  case, 
the  boundary  conditions  are  particularly  difficult  to  predict  since  the 
micro-climatalogical  values  are  quite  unique  for  each  location.  However, 
the  model  is  quite  satisfactory  for  the  evaluation  of  camouflage  materials 
in  that  whatever  boundary  cond  .‘uns  are  used  the  model  will  indicate  the 
response  of  the  candidate  material  relative  to  a particular  leaf 
type. 

From  the  studies  made  with  the  thermal  model  it  was  determined 
that  thermal  infrared  camouflage  is  feasible  with  simple  clothlike  material 
if  the  radiant  properties  can  be  properly  tailored  to  the  background.  For 
example,  it  was  found  that  a material  with  a solar  absorptance  of  0.47  and 
long  wavelength  emittance  of  0.  95  would  emulate  maple  leaves.  It  was 
also  determined  that  this  material  would  not  emulate  leaves  with  larger 
transpiration  rates,  e.  g.  aspen.  These  leaves  are  more  difficult  to 
emulate  but  probably  could  be  satisfactorily  emulated  considering  the 
canopy  of  trees  will  have  large  temperature  variations  due  to  shading. 

In  the  process  of  preparing  the  thermal  model  it  was  noted  that  the 
camouflage  material  must  be  opaque  to  thermal  infrared  if  an  object  is 
to  be  camouflaged.  Any  object  above  the  background  temperature  will 
"shine"  through  a partial  camouflage  net.  This  fact  makes  the  physical 
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construction  of  the  camouflage  net  more  difficult.  A successful  net  must 
have  multiple  layers  of  leaf  sized  camouflage  elements.  These  elements 
must  allow  free  circulation  of  ambient  air  in  order  to  attain  temperatures 
similar  to  the  temperature  of  the  background  leaves.  Furthermore,  the 
material  must  have  good  visual  camouflage  characteristics.  A material 
with  both  satisfactory  visual  and  solar  reflectance  may  be  difficult  to  ob- 
tain due  to  the  overlap  in  these  spectral  regions.  Materials  with  these 
characteristics  have  been  prepared  for  the  near  infrared  but  have  not 
been  reported  for  the  thermal  infrared. 
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QB3MTRL  (S  THE  DETECTABLE  RADIANT  ENERGY  ABOVL  THE  MATERIAL 
SURFACE  IN  THE  3-5  MICRON  RANGE. 

ubbmtrl  is  the  uetectable  radiant  energy  above  the  matepial 
surface  in  the  b-ia  micron  range. 

ub3LEAF  lb  THE  DETECTABLE  RADIANT  tNERGY  ABOVE  THE  LEAF  SORT  ACE 
IN  THE  3-5  M t CRON  RANGE. 

OGbLEAF  IS  THE  DETECTABLE  RADIANT  ENERGY  ABOVE  THE  LEAF  bURFALC 


SURF  ALE 


b-14  MICRON  RANGE. 


; '4*j  U02L 
GvY-'^  L 0?2 


NAMcLIST  / INI/HRS, OT , DECT  / IN2/L AT , DEC , SC , C I , C 4, L ST 
1/IN3/D.KL ,SLP ,AZM  / 1 NA/ ASU  , ASL , A, E , RU , T / IN5/K A, KH, TMAX, T VAR , V 
2/ IN6 /OFS, COND, RG, AG, EG, A SO, RSG 
NR=0 

10  kfcAUl  5,1  NI,tNU  = IOOO) 

READ! 5, IN2 ) 


44 


... , 


rCKlKAN  IV  ^1 


;■* 

■M 


RELEASE  £.0  MAIN 

REA0(5*IN3) 

RE ADt  6 • I N4 ) 

K £ AD  ( 5.1N5) 

READ  (5, IN6) 

NR«NR*1 

GO  TO  ( 2 0 • 30  S « NR 
20  WRITE  (6*107) 

GO  TO  40 
30  wRITE(6*108> 

NK"0 

40  WHITE  (6*1061 
WR I TE ( 6* INI) 
wRITE (6* IN2) 

WRITE (6* IN3) 

WRITE! 6* IN4) 

WRIT  E (6* IN5 ) 

WRITE  (6.1N6) 

C 

C INITIALISE  CONSTANTS 
C 

H*PI 

M*0.0 

ZEN-0.0 

L*0 

1*0 

SINO-SIN(OEC) 

CUSO*COS ( OEC ) 

S I NL  * S I N ( L AT  ) 

CGSL«COS(LAT> 

S I NS" SIN ( SLP ) 

COSS-COS (SLP  ) 

SINA-SIN(AZM) 

COSA-COSIAZM) 

N=o3 • ♦HRS/DT 
WRITEIb.lOOi 
SO  CONTINUE 

H*PI -CST  *0 .262 
CUSH-COS(H) 

S l NH  * S IN  ( H ) 

COSZ*  C'JSL*CUSD*COS  H*-S  INL  *S  I NO 
M*  1 . O/T.O SZ 

IF(M.lT.l.0K.M.GT.6.S)M=0.0 
IMCQSZ.  LT.O.O)COSZ=0.0 
CALL  RAO 
CALL  LEAF 
R T.  S 1 *OA"JR-OC-OE 


FORTRAN  IV  01  RELEASE  2.0 


HA  tN 


00  76 
CO  76 
0077 
007b 
CC7V 
0r>8C 
0081 
CO  82 
0083 
008A 
00  V5 
0086 
008  7 
0088 

0089 

0090 

0091 

0092 
CG93 
009A 
0096 
C0M6 
0'!V7 
009b 
C09S 

0100 
0101 
0102 
010 3 


fN-Tl 
T L“T  L*CV 
oO  KES1-RES2 
70  CONTINUE 
CALL  OETECT 

WRITE  (6,  101)  LST.UATN.GG.USL.USU.QAfQR.QE.QC 
LST«LST*DT/60. 

IF1RL.LT. O.OJGQ  TO  80 
1 «!♦  1 

Xl(l )«LST-(uT/60.) 

X2  ( I ) “TA 
X 3 ( I ) «Tl 
XA(I)«UB3 
X5(I)-QB8 
GU  TO  90 
80  CONTINUE 
i*in 
X6( I l«TL 
X 7 ( l ) «UB3 
Xdt l ) *088 
X9( I )«TGR 

90  IPd.LT.N  )GQ  TO  50 
IHNR.GT.Q)  GO  TO  10 
WRITE  (6  « 10A  ) 

WRITE  (6.103)  (Xl(l)tX2(l)«X3(l)tXb(  1 '.X9( I ) , X A(  l ) , X7 ( I). 
1X5(1).  X8U).I»1»N) 

WRIT  E (6. 105) 

GO  TO  10  

1000  STOP 
END 


✓ 


‘16 

* 

*4 


F C'R  TRAN 

roci 

0002 

0003 
0'JOA 
000*5 
COOb 

0007 

0008 


IV  01 


» I 


coos 

OG  10 


U011 

j 0012 

* CO  10 

S OOlH 

; oo lb 

Mlo 

0017 


RELEASE  2.0  RAO 

SUBROUTINE  RAU 

COMMON/DATAIN/P! .C2. S I G. TL . EPS .G .K2 , CP  * 2 
COMMON  ot.m.lat.dec.h.lst.l.ua.ur.ta.v.tmax.tvar 
COMMON  QB3*U88«Cl»QATM.UU.Ufc»RU»T»SC*C0$TH 

COMMON  OC.HC. ASU.ASL . E . A ,RG  ,EU, A SU ,RSG » AG  .OS T ,QS U .QSL .QDr ,CA . 
COMMON  RH.D.KA.RE.GR.RL. UFS  , CUND. DEL T , TOR 

COMMON  SINO.COSU.SINL.COSL.  SINA.COSA.SINS.CQSS .S  1NH,  COSH,  COS' 2 
REAL  LAT ,l$T.M,KA,K2 
C 
C 
C 

c 

C SUBROUTINE  RAO  PERFORMS  THE  FOLLOWING. (11  ESTIMATES  ThE  SOLAR  ANO 

C UIFFUSE  RADIATION  UN  CLEAR  DAYS  AND  THE  ATMOSPHERIC  AND  GROUND 

C RADIATION. (2 1CACCULATES  THE  RADIANT  ENERGY  ABSORBED  BY  A Lt A h ' 

C OF  ANY  SURFACE  ORIENTATION  AND. (3)  ESTIMATES  THE  AIR  TEMPERATURE 

C HAStO  UN  THE  MAXIMUM  AIR  TEMPERATURE  ANO  THE  MAXIMUM  VARIATION 

C IN  THE  AIR  TEMPERATURE  DURING  THE  SIMULATION  TIME. 

C 

C 

C 

C 

C SULAk  BEAM  RADIATION 
C 

UST-SC/EXP(C<.*M1 

IFIM.LT.  1.0R.M.GT.6.5IQST-0.0 

c 

C DIRECT  SUIAR  TRR AO  I AT  I ON  OF  LEAF  UPPER  SURFACE 

C 

COST  H*S I ND*S INL*CO  SS 
COSTH«COSTH”SlNU*CCSL  + Sl NS* COS  A 
CUSTH*COSO*CGSL*COSS*CCSHtCCSTH 
COST  H«COSD *S  INL*SI  NS*COSA*COSH«-CQSTH 
C OSTH*CUSD*S l NS*S I NA*$ INH*COSTH 
IF(COSTH.LT. 0. 01  CUSTH*0. 0 
QSDT  »CUS  T H*QST 


DIFFUSE  SOLAR  RADI  AT  I ON ( OOF  1 INC  I DE  NT  ON  A SURFACE 

UUF»Cl*UST 

FR AC  * ( 1.0+CUSS1/2. o 

QOFT«u0F*FRAC*UQF*(  l.O-FRAC  )*RSG 

TOTAL  SOLAR  RAOIATIUN  ADSORBED  BY  THE  TOP  SURFACE  OF  THL  LEAF 

OSU«ASUMQSuT+QOFT  1 

AIK  TEMPERATURE  CALCULATIONS 

IML«LST*Pt/l2. 

VAR«5o.29+29.32*CaS<  TME ) *38 .A 8* SI  N( TME 1 -3 . A8*C0S (2. *T ME ) -8  .3 S*S  1 N( 
1 2 • *TME 1 

TA*  TMA  X- ( I TVAR*VAR) /100, 1 

CALCULATING  VAPOR  PRESSURE  OF  820  IN  AIR  AT  TA 
P«(  4.  o*EXP(TA/U.  l*«  1 *AH 


4 


L?t 4?., hisLt-q f I - 


4 


■3 

■ <5 


f Ok TRAN  (V  U1  RELEASE  2.0 


I 00  26 
£ 002  / 
C028 


| on:v 

j Ud  Jd 


- 00  31 

00  32 
003  5 


RAO 

C UATM-  INCIDENT  ATMOSPHERIC  RADIATION 
C 

7AK«TA*2?3. 

OAT M«SIG*(TaK* *♦•!*<• 44+»  08*SUR T( P 1 1 
CALL  UTEMP 
L 

C SOLAR  RADIATION  ABSORBED  BY  THE  LOWER  SIDE  OP  THE  LEAP 
C 

yDPL«OOF*( 1 .0-FRAC1*UDF*FRAC*RSO 
OSL*  ASL  *(  QDFL  ♦•QST*GO  S£*R  SO* F RAC  I 

C 

C TOTAL  RADIANT  ENER GY ( OA ) AbS CRBtD  BY  THE  LEAF 

C 

UA*  A*  l OATM  ♦OG'*-RO*y  AT  M I ♦ Oil  + QSU 

RETURN 

END 


JT*  >V" 


A’ 5 ft  M A^B 


i*  tl 

K «*»*  Cvi* 


I CHTHAN  IV  G1  RELEASE  2.0 


SUBROUTINE  LEAF 

COMMON/DATAIN/PI  ,C2.S  IG.  TL  . EPS.  G.R  2»  CP  . Z 
COMMON  DT*M.LAT*DEC*H»LST*L»OA*URtTA*VfTMAXiTVAR 
COMMON  883,088. Cl.QATM.OG.UE.RU.T. SC.COSTH 

CUMMON  QC.HCiASU. ASL, E .A.RGtEG. ASU.RSG, AC. US T , OSU .0 SL »00F ,CA 
COMMON  KH.DfKA.KEiGR.RL.DF  StCONO.OELT  » TGR  " 

COMMON  SIND.CUSU.SINL.CUSL. SI NA.CO SA , S t NS.CQ SS ,S I NH .COSH.COS Z 
REAL  LAT .1ST .M.KA.K2 


SUDRQUT i NE  LEAF  CALCULATES  THE  HEAT  TRANSFER  TO  THE  LEAF  OR 
MATERIAL  SURFACE  BY  EVAPORATION  AND  CONVECTION  ANU  THE  HEAT  LUST 
BY  THE  LEAF  OUE  TO  RADIATION. 


IF (RL  .LT .0.0 )G0  TO  3 

CALCULATION  OF  KATE  OF  HEAT  TRANSFER  BY  E VAPOR AT  I ON  I UE > 
BOUNDARY  LAYER  EVAP  RESISTANCE 

R8«I60.*0/V»  **0.55 
R8-R2*R0 

TOTAL  EVAP  RESISTANCE 
RT-RL+RB 

DENSITIES  FOR  Q EVAPORATION 


oun 

RHUL«l.O/VG(TLI 

00 1 A 

RHO-i .0/VU(f Al 

00  15 

UE-»RHOL~RH*RHO 

UOLo 

UE*(2«*t>0.*380.Q*0El/RT 

COl  f 

IF  I OE .LT .0.0 ) 06*0.0 

or  ib 

GU  TO  A 

Ool^ 

3 

UE*0  .0 

co?o 

A 

CONTI NUfc 

CALCULATION  OF  THE  RATE  OF  HEAT  TRANSFER  BY  C UNVEC T l UN t cC ) 

TFLM«( TA.TL) /2.Q 
ANU*B.032«-A,8622E-2*TFIM 
ANU»ANU  + t..06E-5*TPLM*TFlM 
RE  «V*0/ANU 

oK»Ci*U**3*ABS(  TL-i  A)  / ( ANU*ANU*TFLM) 

IF ICR/RE**2J0.LT .1.1  GO  TO  l 

CALCULATION  OF  FREE  CONVECTION  COtF.(HC). 

HC l * ( 0. A V 7*KA*GP  **0. 251/0 
HL2»HCl/2.  * 

HC*HC1+HC2 
OC=HC*(  TL-TAI 
GO  TO  2 
l CONTINUE 


I- OK  TRAN  IV  (,l  RELEASE  2.0 


DETECT 


SUBROUT  1 NE  DETECT 

DIMENSION  B150),FV(50),WAVE150> 

CUMMON/O AT  AI N/Pl , C2 , S IG , TL , EPS.G ,K2 , CP , 2 
COMMON  0T,M, LAT . OEC. H, LST ,L  ,QA,  OR,  T A,  V.  TMA X,  TVAR 
COMMON  0B3,0B8,C1 , QA T V , OG, UE , RU , T , S C , CQS1  H 

COMMON  UC.HC ,ASU,ASL ,E ,A,RG  ,EG  , A So  , R SG  . AG , CST  , OSU  .QSL  , OOF  , CA 
COMMON  RH.D.KA, RE , GR, XL. OFS . CUNO, OEL T , TGR 

COMMON  SiNO, COSO ,$INL,COSL,SlNA,COSA,S INS,  COSS  ,S  INH,C0SH,C0S2 
REAL  LAT,LST,M,KA,K2 


THIS  SUBROUTINE  CALCULATES  THE  RAUlANT  ENERGY  WHICH  WILL  BE 
DETECTED  ABOVE  THE  RADIATING  SURFACE  IN  THE  3-6  AND  THE  8- 1A 
MICRON  RANGE 


WAVE  11) = 3. 
wAVcI  2)  = 5. 

WAVE  C 3 ) = a . 

WAVE ( A)-l A. 

T LK=TL*273. 

UO  3 1=1, A 

Bin  =C2/(WAVE(l)*TLK) 

IF(B(I).LI.2.)GJ  TO  A 
FVII )=0.0 
00  5 N = 1 O 

FVl  I)=Fvi  1)+.15*(((EXP(-N*8(l))) /N**A. )*l ( IN* B ( I ) *3 . ) *N*b ( I ) *b . ) 

1 *N*B( I ) *6 . ) ) 

5 CONTINUE 
GO  TO  3 

A FVI  IJ-l.-U.  I S*B(U**3.)*t  .333-dll  )/6.-Mb(  I J ) **2  . /bO.  - ( B ( I ) I ** A. /6 
1 AO,  + ( B ( 1 ) )**6.  n 72  16  0. -(<31  I >)**B./  13305600.1) 

3 CONTINUE 

Fl«FV(2)-FV(l) 

F 2 = F VI  A)-F  V(  3) 

F3= .016 
F A = • 0 0 2 

Urt3=( u*F  1*  SI G* l T LK**  A. ) ) ♦( F 3*RU*i . 9A*C0ST H) ♦ I F 3 *T +RSG* 

1 1 .9A*C0S  L ) 

Odb=(  E»F2*'bI  G*  (TLN**A  . ) ) » I FA*RU  *1 . 9 A*COST  H ) ♦ I F A*T  *RSG* 

11. 9A*C0SZ ) 

RETURN 

END 


w v'.‘4 


i£m. Vf  y 


•sa?j»*5a»ss; . • 5W--JLiy4»  ..1; l 


M ft>v 

tif,. 

m w- 

|v  r. 

P fe 

i m 

T;  %i 

, 

g\  IS; 

w p 

I oh  Than 

1# 

|a  ■ 0 

ooo  i 

pi: 

CO  02 

00  0 2 

i-'l'' 

' ODD  A , 

|-  % 

ooos 

»;  S- 

OOOo 

% 

000/ 

1 f 

00  * J 

ooo'v 

RELEASE  2,0'' 


GTEMP 


SUBRGUT InE  GTEMp 
DIMENSION  TG  (200 ) 

C.OMMON/JATAIN/PI ,C2, SIG, TL , EPS ,G , K2 , C P ,Z 
COMMON  OT.M.LA.T, DEC, H.LST.L.QA, UR,  TA.V.TMAX.TVAR 
COMMON  OB3,U88,C1,UATM,OGjUE,RU,T,,SC,COSTH 

COMMON  UC.HC  • ASU.A  SL.E  »A  ,RC  »EG  » ASG  .RSG.AG.GST  » USU  » OSL  » QDf,  Ca 
COMMON  RH.O.KA, RE,  oR.RL, OPS, CbNU,L)L-l.T,,TGR 

COMMON  SI  NO. COSO, S'  INL.COSL.StNA.CuSA.S  INS.COSS.S  INH.COSH.COSZ 
REAL  LAT ,L  ST.M .KA.K2 


SUBROUTINE  GTEMP  CALCULATES  THE  GROUND  SURFACE  TEMPERATURE 
BY  AN  EXPLlCT  FINITE  DIFFERENCE  Me  THUD. 


N IS  THE  NUMBER  CJF  NODES  UStD  IN  CALCULATING  THE  TEMPERATURE 
GRAOIENT  IN  THE  SOIL. 

UELTX  IS  THE  DISTANCE  BETWEEN  NODES. IN  CM. 


oulo 

N*20 

coil 

’ TI  ME  -0.0 

oo  i 

DtLTX=.SOKT<  2.  +OP  $*OELT  ) 

OCI  J 

CSI-QST+C0S2 

001*. 

CsI  ^OS  I «•  CDF  . 

00  1 S 

IF(C.GT.O)  GO  TD  2 

0(1  lo 

‘ 

JJ*N-l 

L01  7 

DO  3 i *1.  • N 

JulO 

3 

TGI  I ) = TMAX-TVAR/2»  0. 

oolo 

2 

CONTINUE 

00  20 

TGI*  TGI l X 

(i021 

C5*0F  S*U6.'.  T / ( UEL  TX  *0E  L t X 1 

CG?2 

C6--CGNu/0ELTX 

OO  2 3 

* i 

TGK*TGl  1X2  7 3.  0 

co?-; 

’ UG*EG*S IG*TGK**A 

0u2  b 

TFLM*(TGl«-TAI*0.5 

002c 

ANU*  d.  08  2*  A. 8o22E-2“TF  LM 

0027 

ANU* ANU+8 . 06 t~ 5* TFLM* TFLM 

C02B 

RHC1A*0 .3537  (TA+2/3  .0) 

Ou2s 

PR*  1. 0/(  CP*RHUA#ANU/*A  )**0.  66 

00  30 

CF*1  .0/  < 2.5*ALOG(  lfeO.O/t  1*5.0)**2 

Cnil 

H(j*PR*CF*RHOA*CP*V 

00.0  2 

UUAn=wS  I *ASG+wATM<‘AO-UG 

003.1 

TG  (11  = I HG*T  A*-C6*Tb(  2 ) +3UAN  1 /IMGfCfe  i 

00  2 A 

00  i 1 *2 , J J 

JO  30 

l 

TGI  I } = TGU  X-CS*ITG(  l + U-2.0*TGl  I I^TGI 

Co  3<j 

IF  (L.C.T  .0  ) GU  TG  A 

00  j / 

IFlAbSl  TGl-TGl  m.LT..30B)  GO  TO  S 

or 

GO  TO  2 

C03S 

A 

CONT  I NUE 

''OAc 

T i ME  * TI ME ♦UF  L T 

JOA  1 

IKT1ME.LT  .ot  1 60  TO  2 

U0A2 

L = L*1 

OJA  3 

Ti.R*.TG(  1) 

JOAA 

RETURN 

Lv-Ab 

c ND 

■-  ‘iitasSsWii&sA'  ‘aai  jjjjg ift  Ysv  t arils'  it- 


BLK  DATA 


✓ 


S3 


1=  OP. THAW  IV  U1  RELEASE  2.0 


I cao  i 
0u02 
0003 
000  A 
OuOS 


Cu06 


t*'  OCK  DATA 

LOKMON/OATAIN/Pt ,C2, SIG, TL  , EPS.G ,K2 ,C P . Z 
REAL  LAT  ,LST ,M ,KA,K2 

DATA  PI  ,C2  .SIG  ,TL,EPS  ,G/3.i.Mi6,  143 BS.  , . 822  E- ID  ,1  0 .,  .005.  .35  3 EOT/ 
DATA  K2,CP,Z/.062,.24,5./ 

C 

c 

c 

r.  SUBROUTINE  BLOCK  DATA  ALLOWS  VARIABLES  IN  COMMON  TO  INITIALIZED 

C IN  A DATA  STATEMENT.  T HE  FOWWOW  ING  DATA  IS  INITIALIZED  AS  BLOCK 

C DATA. 

C C2  IS  A CONSTANT  USED  IN  PLANCKS  SPEC  TP.AL  ENERGY  DISTRIBUTION 

C .IN  NIC  RONS -UcGRC  ES  K 

C SIG  IS  THE  STEFAN-BOLTZMAN  CONSTANT, IN  CAL/CM**2 -K**V 

C TL  IS  THE  ARBITRARY  INITIAL  TEMPERATURE  OF  THE  LEAF  OR  MATERIAL 

C SURFACE. 

C.  EPS  IS  THE  PRECISION  WITH  WHICH  THE  TEMPERATURE  OF  THE  LEAF  OR 

C MATERIAL  SURFACE  IS  CALCULATED.  NOTE,  IF  EPS  IS  SET  100  SMALL 

C A DIVIDE  CHECK  WILL  OCCUR  IN  THE  MAIN  PROGRAM  DU E TO  THE  SMALL 

C CHANGE  OF  RE  SI  AND  R6S2 

C G IS  THE  GRAVITATIONAL  CONSTANT, IN  CM/MIN**2 

C k2  IS  A PROPORTIONALITY  CONSTANT. 

C CP  IS  THE  SPECIFIC  HEAT  OF  AIR, IN  CAL/GM-UEGREE  C. 

C L IS  THE  GROUND  COVER  V EG  AT  ION  HIEGHT,IN  CM. 

C 

C 

c 

END 


APPENDIX  B 


Evaluation  of  Various  Surface  Coatings  1 

■ i 

on  Camouflage  Material  Temperatures  | 


Introduction 


As  shown  in  Table  1,  camouflage  material  temperatures  are  several 
degrees  higher  than  foliage  when  exposed  to  solar  radiation  and  slightly 
lower  during  the  evening  hours  when  no  solar  insolation  is  present.  This 
information  leads  us  to  the  conclusion  that  the  radiative  characteristics 
of  the  camouflage  material  must  be  altered  if  it  is  to  emulate  foliated 
backgrounds.  To  accomplish  this  end,  two  basic  approaches  were  con- 
sidered: (1)  change  the  radiative  characteristics  of  the  top  surface  of 
the  material  by  using  a clear  spray  coating  or  a transparent  acetate  cover 
as  previously  discussed  and;  (2)  alter  the  solar  and  I.  R.  energy  absorbed 
and  emitted  from  the  lower  side  of  the  material  in  order  to  control  its 
temperature. 

Two  types  of  surface  configurations  were  used  to  evaluate  the  second 
approach.  The  first  consisted  of  bonding  a sheet  of  aluminum  foil  to  the 
lower  surface  of  the  camouflage  material.  The  foil  has  the  effect  of 
lowering  material  temperatures  when  large  amounts  of  solar  radiation 
are  reflected  from  the  ground  and  has  little  or  no  effect  when  no  reflected 
solar  radiation  is  present.  This  is  due,  in  part,  to  the  fact  that  the  foil 
effectively  eliminates  the  absorbed  :short  wavelength  radiation  which  is 
reflected  from  the  ground,  and  eliminates  the  long  wavelength  exchange 
between  the  lower  surface  and  the  ground.  In  addition,  another  piece  of 
camouflage  material  was  coated  on  the  lower  surface  with  white  lacquer. 
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The  lacquer  has  the  effect  of  reducing  the  solar  radiation  absorbed  on  the 
lower  side  of  the  material  while  the  long  wavelength  emittance  remains 
unchanged.  Thus  the  white  paint  effectively  lowers  the  temperature  during 
periods  when  solar  radiation  is  present  and  has  little  effect  at  night. 


Experimental  Proceuv. rc . Results  and  Conclusions 


Temp  .ratures  of  the  plain  camouflage  material  and  four  variations 
thereof  were  made  and  compared  to  a Botanical  Wonder  plant  (Fatsia 
Japonica).  The  measurements  were  made  with  Barnes  PRT-5  and 
PRT-10  radiometers  and  were  carried  out  with  varying  atmospheric 
conditions  in  order  to  properly  evaluate  the  effect  of  the  coatings.  The 
results  of  these  measurements  are  presented  in  Tables  B-l  and  B-2.  In 
addition  the  air  temperature,  plant  temperature,  plain  camouflage 
material  temperature  and  one  variation  of  the  plain  material  temperature 
were  plotted  using  data  from  the  PRT-5  and  are  shown  in  Figures  B-l 
through  B-4. 

No  concrete  conclusions  can  be  drawn  from  the  preliminary  data 
obtained  thus  far;  however,  certain  trends  are  evident.  The  white  lac- 
quer coating  lowered  the  material  temperature  below  that  of  the  plant  and 


plain  material  during  daylight  and  evening  hours  (Figure  B-3).  Both 


transparent  coatings  effectively  lowered  the  material  temperature  during 
daylight  hours,  however  the  spray  coating  increased  the  temperature  at 
night,  while  the  acetate  cover  tended  to  lower  the  apparent  temperature 
at  night.  Unfortunately,  the  acetate  cover  reflects  large  amounts  of  short 
wavelength  radiation  thus  producing  glare  (Figures  B-l  and  B-4). 


! -I 


iriTiOOOOinLriinirooinoLn 


UM 
• ^ w 

W * s 
a ^ o 
p > 


oo>riooooooou‘>tjr»iriO!jnoooiriOL^oooc. 

pH  —*4  O0  pH  —4  —4  fSJ  00  ^4  ^H  —4  —4  pH  i—4  —4  i— 4 pH 


r 

> 

0 

'J 

c 

O 

lTi 

0 

0 

0 

0 

O 

O 

LT 

in 

m 

m 

0 

© 

O 

O 

O 

0 

0 

in 

0 

U 

cq 

r- 

r\3 

M 

sd 

iA 

00 

M 

CO 

m 

00 

00 

M 

OO 

N 

N« 

ni 

(N) 

CS3 

00 

t 

»“4 

•— * 

<NJ 

(N 

ro 

00 

*-♦ 

f— * 

p«H 

*W 

^■4 

«-*4 

p—4 

«“4 

• <*4 
0 
U4 

c 

0 

minooooir.  ooooooooooooriiririOLno 

NNOtO'OOO^lM'iiy'N^PJMOOHHKir(imNrf,i4j 

■— < PsJ  4*-H  »*H  (NJ  (\J  fSj  m-i  r-^  p— l M *— 4 —4  *"*H  r«^  P-H  *■»* 


oooooooooooouooir  oomo>riLr'OOoin 

^“4  nJ  nJ  ^4  OsJ  rvj  f\J  nj  "H  r— 4 »“4  »»«4  «"»4  »— 4 


.3 

‘;j; 

IV 

y 

■ij1 

. ;» 

4* 

» ..' 

V 

r S " i r » 1 1 

! j] 

41 

B&BH 

* 

f 

ss® 

..  rg&Sfi 

- .»•-•• 

« S ° 

s rs 


cJ  U 
u u r 

•-i  fl) 

< 2-.S 

f! 

4) 

H 


OOOOOOUOLfiOOOOOOunCOOOinOOOOO 

Oirt'hN^NM'tff''£NNN'tMHMW(fii\jH-4-i~4 

^4  «~4  P"-4  psj  fSJ  ^sj  p ■ A **«4  "“4  P~‘ 


uoif'OOOOOiriiPOOi.riiriAir'iriOOinc  o u*.  m in  © 

■+^■10  040^40000- i o c © © o © © © m in  r*  «+  -t"  ■<* 


oooo’^inor'ooo'o— i oj  w ^ in 

hhhhhhh-'MNNNMN 


I I I I I I I 


oooooooooooooooooooooooooooo 

4 ««•«»«•  ••«•••««  «•  * ?•  * ' % • •» 

NO\oNf~,fooo'^t''^,'i,Lnooosr-r-oor~'Or-i.r'Aj'rt,'!f'^'iriTj*pri 


OOOOOOOOOOOOmOOOOOOOOOOOO 
't  O O'  "t  Tf  od  O-  00*  o’  O O O O O'  O'  oc  oo‘  o’  — <M*  — * t\j  M*  tA 

-N-tNNNCl-iNH-iwrtH  I i i i i i 


coocoooooooooooooooooooooooo 

r0f0h~c^O'0('i—»00'Oif'OC'0'00000'r^r^r~s0P0p0f^f^i'j'('jN 

- <-<-rfNHNNNN-‘M  I I I I I I.  I 


OOOOOOOOOOOOOOOOOOOOOOOOOOOC 
TCUlt'WO'flN- i00in'0O0'0'03000't^^-r'r}<(Mr^f^'y»l<r''vj 

I I I I | I I 


ooooooooooooooounoooooooooooo 

m^o>oooa'HO'f^o>o'NOOO-oooc'NKimi<u»i‘<tw 

~ <l— ^ ^^^hNnnnhhh  >-«  m — I — I -J  I i | | i i | 


inmin©©o©ir'uiu'iri©om»rmininooi/'>cooiftiri»r>o 

't'f'fiAO^'O'O'O'O'Or'- •00000©0'T'K'U'inf'rf"3,-t 


— i-j-j  — — ‘ — i— 'MMMNrONNNN 


DATA  POINT 


DATA  POINT 


DATA  POINT 

Figure  B-4 


63 


* " • .. 
fc 

S? 

*3  ■ - . - ' . 

V; 

The  foil  cover  produced  the  most  promising  trends.  It  effectively 
p-  lowered  temperatures  during  the  mid-morning,  and  mid-afternoon  hours 

f;  and  slightly  raised  temperatures  during  the  evening  hours.  In  addition, 

it  followed  the  plant  temperature  more  closely  than  the  other  variations 
- (Figure  B-2).  This  configuration  holds  the  ad  ted  advantage  that  most  of 

the  radiation  emitted  from  a hot  object  placed  under  the  material  would 
be  reflected  off  the  lower  side. 

The  measurements  made  thus  far  are  only  preliminary  and  were 
designed  to  establish  various  trends  by  altering  certain  radiative  pro- 
perties. It  is  not  known  if  the  materials  tested  could  be  used  under  field 
conditions,  however  it  is  evident  that  progress  can  be  made  towards 
emulating  foliated  backgrounds  with  continued  research  in  this  area. 
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